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Abstract 

We introduce a simple mean- field lattice model to describe the behavior of nematic elastomers. 
This model combines the Maier-Saupe-Zwanzig approach to liquid crystals and an extension to 
lattice systems of the Warner- Terentjev theory of elasticity, with the addition of quenched random 
fields. We use standard techniques of statistical mechanics to obtain analytic solutions for the full 
range of parameters. Among other results, we show the existence of a stress-strain coexistence 
curve below a freezing temperature, analogous to the P-V diagram of a simple fluid, with the 
disorder strength playing the role of temperature. Below a critical value of disorder, the tie lines in 
this diagram resemble the experimental stress-strain plateau, and may be interpreted as signatures 
of the characteristic polydomain-monodomain transition. Also, in the monodomain case, we show 
that random-fields may soften the first-order transition between nematic and isotropic phases, 
provided the samples are formed in the nematic state. 
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I. INTRODUCTION 



Liquid-crystalline molecules imbedded in a polymer network give rise to the novel class 
of elastomer systems, with coupled rubber elasticity and orientational order, and rather 
unusual properties l|, |2| . Nematic elastomers (NEs) may undergo a distortion in response 
to an alignment of the nematic units as the sample is cooled below the nematic-isotropic 
transition temperature Tni. Reciprocally, the application of an external stress may give rise 
to nematic ordering of an initially disordered sample. There is a number of suggestions 
of applications of these new soft-matter materials, ranging from uses in optics (in bifocal 



lenses, for instance) to applications as thermo-mechanical devices 

Compared with conventional liquid crystals, NEs present a peculiar transition from the 
nematic to the isotropic states. In usual nematics, according to the Landau-de Gennes 
theory, sjTiimetry requirements lead to a first-order transition, with a jump of the nematic 
order parameter, at a transition temperature Tni. A different scenario is observed in the NEs. 
Instead of a discontinuity there is a continuous but quick variation of the order parameter 
during the transition. The microscopic mechanism behind this non-trivial behavior has 
been discussed by many authors js-?]. In particular, continuous three-dimensional coarse- 
grained theories |(|, and numerical simulations of microscopic models {s], indicate that 
quenched random-field interactions may smooth out the characteristic first-order transition 
in these systems. In NEs these random fields are supposed to originate from internal stresses 
produced by the network cross-links jsl, More recently, slightly different random-field 
interactions have been ^ggested by Lu et al. to describe network heterogeneity in randomly 
crosslinked materials 



Random fields may also be relevant to describe the interesting transition from states of 



polydomain to a monodomain in NEs ll|. If the system is cooled below T^i, NEs are known 
to display a stable structure, the so-called "Schlieren texture" , characterized by sets of frozen 
randomly-oriented domains of mesogen units. This polydomain state may be turned into a 
monodomain state by stretching the sample under uniaxial tension. For a certain range of 
low temperatures, experiments show an unusual stress-strain curve with three characteristic 
regions. In the first region, for small strain, the system is in the polydomain state, the 
sample is opaque, and stress increases linearly with strain, according to Hooke's law. In the 
second region, the stress is constant for a range of intermediate values of strain, and there 



is a dramatic increase of the nematic order parameter. This plateau is then followed by an 
increase of stress at larger strains, with the system in the optically transparent monodomain 
12-|l4|. authors have suggested that cross-linking conditions are essential to 



state 



explain the stable polydomain state 



In particular, random-fields and the neoclassical 



theory of elasticity have been used by Fridrikh and Terentjev to obtain a good fitting of 
some experimental findings 11 1. 

In this article we propose a lattice statistical model for nematic elastomers, at the mean- 
field level, which can be analytically solved in the presence of stress and random fields. 
This model has the advantage of providing a simple way to investigate a wide range of 
parameters, gives a unified view of the critical behavior of NEs, in general agreement with 
previous theoretical jo]. 111, numerical jsj, and experimental findings l|, and suggests the 
occurrence of some additional phenomena. In Section II, we define the basic lattice model. 
The global free energy is obtained in Section III. Some specific calculations, including stress- 
strain curves, comparisons with the literature, and a few new predictions are presented in 
Section IV. The main conclusions are given in the last Section. 



II. LATTICE MODEL FOR A NEMATIC ELASTOMER 



The Maier-Saupe (MS) model, which is known to provide a good description of the 



16|, is the liquid-crystalline analogue of the Curie- Weiss 



isotropic-nematic transition 

model of ferromagnetism [l^. The idea consists in the extension of the finite-range interac- 
tions of realistic systems to infinite-range interactions in order to construct a simpler model, 
which turns out to be exactly soluble. In the MS model, the basic elements are the molecular 
unit vectors nj, representing a set of mesogens, which interact via the quadrupole energy 

J, _ A 

where A > is an energy constant. 



l<i<j<N fi,u=x,y,z 



(1) 



S: 



flU 



-(3nr<-5n 



(2) 



and \ni\ = 1 for i = 1, 2, ■ ■ ■ , A^. The calculations may be further simplified if we adopt a 
suggestion of Zwanzig and restrict the number of allowed orientations of rij to the six values 
along the Cartesian axes, 

G {(±1,0,0), (0, ±1,0), (0,0, ±1)}. (3) 



Generalizations of this model, which we call the Maier-Saupe-Zwanzig (MSZ) model, have 
been applied to the study of biaxial and uniaxial ordering in rod-plate mixtures of liq- 
uid crystals 18|, |l9|. In agreement with predictions of Landau-de Gennes theory 16|, the 



isotropic-nematic transition is found to be first-order, with a discontinuity in the nematic 
order parameter. In the nematic elastomer case, Xing et al. have considered similar Maier- 
Saupe type of interactions adjunct to a microscopic model to make contact with a macro- 
scopic Landau theory [2^. 

character of nematic elastomers, a number of 

approaches have been considered 



To model the mechanical and orientational c: 
arise from an entropic contribution Thus we consider the canonical partition function, 



23| . We shall assume that elastic properties 



Z = J] l]({n,},A)exp(-/3E, 

{rii} 



MS J 



(4) 



where /3 = l/fc^T, Ey[^ is the interaction energy of the MSZ model given by Eq. ([T]), and 

the sum is over the configurations ([3]) of the microscopic nematic directors nj. The entropic 

;erm VL depends on the nematic orientations and on a global lattice distortion tensor A 
1 

1[. If we consider a uniform strain along the direction of a unit vector m, the distortion 
components may be written as 



(5) 



According to an extension of the neoclassical theory of elasticity |2|] for lattice Hamiltonian 
systems, proposed by Selinger and Ratna [5i], we write the degeneracy as 

n = exp (-/3Fei) , (6) 



with the elastic free energy 



^ei = fETr(^o.-A^.Zri.A) 



(7) 



where is the linear shear modulus, li is a local shape tensor, and Zo,i is the local shape 
tensor at the time of the cross-linking. The components of the shape tensors are obtained 
from the equation 

= ^I^^af} + (l^^ - ^I^) ni^a'ni,l3, (8) 

where l± and l\\ are the effective step lengths of the nematic polymers in the perpendicular 
and parallel directions with respect to the nematic vectors. If the cross-linked network is 



formed with the sample in a totally disordered isotropic state, we assume that the shape 
tensor l^j is given by an isotropic average of 



--( 



2/ / + /„ M (5q,/3- 



(9) 



Thus we obtain 



N 

E 

i=l 



(A^ + 2A-)-5(A^-A-^)(^(m.n.)^-l 



where we have introduced the tensor 



^ ' i=l u.u 



(10) 



and the parameter 



1 o;-l 



6 



217' + ' 



0<6 <1, 



(11) 



(12) 



with 5 = in the isotropic case, and 5 = 1 in the limit of largest anisotropy. The first term 
on the right-hand side of equation ffTOl) is the classical rubber free energy. 



(13) 



Finally, the partition function may be written 



Z = ^exp(-/5£'. 

{rii} 



eS 



(14) 



with the effective energy 
A 



E. 



cS 



N 



(15) 



l<i<j<n fi,iy=x,y,z 



1=1 fj,,u=x,y,z 



where we have defined 



(16) 



Given the entropic origin of the elastic contribution to the free energy, the linear shear 
modulus should depend linearly on temperature. 



(17) 
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2l| 



where Ug is the number of strands in the polymer chain per unit volume [l|. 

We further assume that the microscopic nematic directors are subjected to random fields, 
which originate from the distribution of local anisotropy axes generated at the time of cross- 
linking. These interactions may be represented by the energy term 8|, |9 

N 

i=l fi,i'=x,y,z 

with 

Hr = \ {^KK - , (19) 

where 7 is an energy parameter, and is the fi-th component of the unit vector hi. We as- 
sume that {hi} is a. set of independent and identically distributed quenched random variables 
with probability distribution 

. . f c/2, for h = (0,0, ±1), , , 

P{h) = < V , , ;> ^20) 

[ (l-c)/4, for ^ = (±1,0,0), (0,±1,0), 
with < c < 1. Let us choose the symmetry axis of the mesogen units along the z direction. 
Then the parameter c is related to the degree of anisotropy of the mesogens attached to the 
network cross-links, so that c = 1/3 for samples formed in the isotropic state, and c > 1/3 
for samples formed in the nematic state. Note that we have assumed the shape tensor 1^1 be 
fixed as an isotropic average of and chosen to include the network- heterogeneity history 
in the random-field interaction only. 

III. CALCULATION OF THE FREE ENERGY 

Effects of a fixed external stress a will be taken into account by means of the partition 
function 

dX e^-^" Yl (^eff + E,i)] , (21) 

for a given a configuration of random fields. The free energy density is given by 

/ = lim ^ In F = /rub - - ^ lim ^ In V exp (-/3E) , (22) 

{ni] 



where / should be a minimum with respect to A and 



TV \ ^ AT N 

^ = E I E - 9^ E E ^-^r - i E E ^r^r, (23) 

ii,v \ i=l / i=l iJL,v fi,i/ i=l 



where we have discarded terms of order smaller than N. 

Using a set of standard Gaussian integral transformations 
tions between different particles, 

NI3A 



17] we decouple the interac- 



J2 exp i-PE) = I [dQ] exp ( J] Qj, 



(24) 



where [dQ] = Yl^u V PAN/2-KdQ^j,y. Performing the sum over the orientations n of a single 
particle we obtain 

N(5A 



/ = /rub + |-^A-iln2-rMim iln j [^/Q] exp (^E^m^ + ^^Q^ 



N 



exp<^^ 



4 



TrHi + ln J]e^(/i, 



where 

e^{hi) = exp 
Invoking the law of large numbers, we have 



(3 ^—Q,, + -fHi 



(25) 



(26) 



JL^. ^ E {-T ^^^^ + ^-(^^)) } = "T ^).+ (in (e e.(^) 



(27) 



where denotes the expectation value with respect to the random-field variables, from 
which we see that the free energy is self-averaging. Carrying out the integration using 
Laplace's method we arrive at 



B 1 
/ = /rub + — - c^A - - In 2 - max L{Q 



(28) 



where Qnu maximizes the functional 



PA 



(TrQ2 + TrQ) + (ln (^E^^j^ 



(29) 



The condition for L to be stationary with respect to Q^^ leads to the equations of state 
for the order parameters 



Q 



1 



3 ' ^ 



a^a / h 
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(30) 
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Notice that TrQ = 0. The condition for the free energy to be stationary with respect to A 
leads to the equation of state for the distortion, 



(31) 



where we have used the result (ISUil . Using these equations of state we may rewrite the 
free-energy density as 



/ = /rub + ^ - - iln2 + 4 Tr Q2 - i (in ( ^ 



/3 




(32) 



To make a closer contact with experiments on liquid crystals, we use the standard diagonal 
parametric form of the traceless matrix Q appropriate for the nematic ordering along the z 
direction, 





2 

S 



\ 








(33) 



The nematic order parameters S and ri characterize the isotropic phase (S* = = 0), the 
uniaxial phase (5* 7^ O,-?] = 0) and the biaxial phase (5* 7^ O,-?] 7^ 0). From the equation of 
state fl30D we find 




1 

2' 



(34) 



and 



3 / C-r 



(35) 



The simplest version of the MSZ model defined by Eqs. ([UlS]) lacks the proper symmetry 
to describe a stable nematic biaxial phase. To check if this behavior is robust after includ- 
ing elasticity and random fields, we have considered the case where the strain direction is 
perpendicular to the chosen axis of symmetry of the order parameter. In fact, there has 
been a lot of debate in the literature about the soft (or semisoft) response of NEs being 

nnn 

related to the onset of biaxial behavior [Jj, l24l-l26| . While the description of such phenomena 
is beyond the scope of this paper, we emphasize that a numerical inspection of the equations 
of state have led us to conclude that biaxial symmetry remains absent for a large range of 



8 



parameters. Hence we shall take rj = and consider the strain direction parallel to the z 
axis, m = (0,0, 1). This is a reasonable assumption since the coupling between elastic and 
orientational degrees of freedom provides an easy axis of symmetry for this system. Thus 
the free-energy density fl521) becomes 

1 1 Q/l 1/ / \\ 

(36) 



where 



e^,{h) = exp |/3 
e^{h) = exp <^ /3 



\aS + ^7 {Shi - 1) + 5 



(37) 
(38) 

and the equation of state for distortion (I3T1) takes the form 

^-F^^K''^^)^- '''' 

The results obtained thus far are valid for arbitrary random-field distributions. Hence- 
forth we limit ourselves to the discrete distribution (!20l) . In this case the free-energy density 
( l36l) is given by 



/ = /rub + ^fi - - i In 2 + ^ASiS + 1) 

{cln [2 + e/5(9^^?/4+97/8+i?)] + (1 _ c) In [l + e^'^^/^ + 6^^(9-4^/4+^)] | 



(40) 



and the equation of state for the order parameter flMj) becomes 



^"2 



/3(9AS/4+97/8+B) 



+ 



c)e 



/3(9AS'/4+B) 



2 + g/3(9AS/4+97/8+B) ' ]_ _^ g9/37/8 _^ g/3(9AS/4+B) 



(41) 



IV. THERMODYNAMIC RESULTS 



Without loss of generality, we assume Ug = 1 and 6 = 0.5, since other choices lead to 
qualitatively similar results. Let us write down the equations in terms of which we performed 
numerical calculations. We express the energy in units of A and the temperature in units of 
ks/A. The free-energy density fHT]) is given explicitly as 



/ = ^S{S + 1)+t\ 1{5X^-7X-^) -ln2-cln 



2 + e 



9(2S'+7)/8T+3(A2-A-1)/8 



c) In 



1 + e^'^/®'^ + e9^/4^+3(^'-^"')/8 



(42) 



The equation of state for the order parameter (HTi) takes the form 

^9(2S+7)/8T+3(A2-A-1)/8 Q _ ^\g95/4T+3(A2-A-i)/8 



S 



+ 



2 + e^(2S+7)/8T+3(A2-A-i)/8 ' ]_ _^ g97/8T _^ g9S/4T+3(A2-A-i)/8 

and the equation of state for the distortion (139!) gives 



1 

2' 



(43) 



A 



1 

A2 



a 
f 



2A 



1 

A2 



S. 



(44) 



In the absence of apphed stress (a = 0), the equation of state (l44l) can easily be solved 
for the distortion with the result 



A 



1 + 5/4 



(45) 



,1 - 

The distortion increases monotonically with S, for < S* < 1. In addition, 5 = implies 
A = 1 , indicating that the strain is driven by the orientational ordering only. For an arbitrary 
applied stress (a > 0) the equation of state (jUj) is a cubic equation in A with only one real 
and positive root given by the formula 



A 



2a 



3(2 - S)T 



1 + 2 cosh 



cosh"^ ( 1 + 



27(2 - Sf{A + S)T' 
32^3 



(46) 



which is a monotonically increasing function of S. 

In our numerical calculations we solved Eq. (I43l) for the order parameter S using for A 
the result given by Eq. (H6l) . 



A. Results in the absence of disorder (7 = 0) 

In Fig. [T] we plot the order parameter S and the free-energy density / (upper curve) 
as a function of temperature for applied stress a = 0.02 in the temperature range where 
a first-order transition takes place. The low temperature stable solution abc with larger 
order parameter S will be called nematic solution, whereas the high temperature solution 
dbe with smaller order parameter will be referred to as isotropic solution, even though S 
no longer vanishes as in the case of zero applied stress. The branch cd corresponds to an 
unstable solution. In the narrow temperature interval < T < both the nematic and 
the isotropic solutions are stable, and we have to choose the one with smaller free energy. 
The first-order transition between nematic and isotropic phases occurs at the temperature 
Tni = Tb where the free energies become equal. We observe that, in agreement with classical 

10 
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0.88 



-0.28 



-0.3 f 



-0.32 



0.96 



FIG. 1. Order parameter S and free-energy density / (upper curve) as a function of temperature 
for 7 = and applied stress a = 0.02 



elasticity theory, the free energy changes hnearly with temperature in the isotropic phase 



In Fig. |2]we plot (a) the nematic order parameter S and (b) the distortion factor A as a 
function of temperature for several values of the applied stress a. As expected from Eq. (H6l) . 
the graph of A follows closely that of S. For small applied stress a the system undergoes a 
first-order transition with a gap between the nematic and isotropic solutions. As the aligning 
stress a increases, the gap decreases until the critical point is reached, beyond which there is 



27| 



no phase transition. This behavior has been predicted by de Gennes in the mid-seventies 
before nematic elastomers were proven to be chemically feasible. Experimentally, however, 
no first-order transition is observed down to the limit of zero applied stress. The jump in the 
first-order transition is smoothed out, being replaced by a continuous but quick variation 
of the order parameter. This fact has been interpreted theoretically as being due to the 
anisotropic distribution of random fields js, 6|. 
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FIG. 2. Nematic order parameter (a), and distortion factor (b) as a function of temperature, for 
7 = and various applied stresses a. 

B. Effects of disorder (7 > 0) 

Let us examine how the random fields affect the nematic-isotropic transition. In Fig. 
we show the nematic order parameter as a function of temperature for cr = 0, 7 = 0.2 and 
several values of c. We observe that the gap between the nematic and isotropic solutions 
decreases as c increases from c=l/3toc = 0.422, disappearing above this value of c. In 
Fig. [3]d the nematic order parameter is shown as a function of temperature for c = 0.44 and 
several values of 7. Again, the gap between the nematic and isotropic solutions decreases as 
7 increases from 7 = to 7 = 0.263, and the transition disappears above this value of 7. 




1.05 




FIG. 3. Nematic order parameter as a function of temperature for a = 0. (a) 7 = 0.2 and various 
values of c. (b) c = 0.4 and various values of 7. 



These results show that anisotropic (c > 1/3) distribution of the random fields of sufficient 
strength (7 > 0) is necessary to smooth out the isotropic-nematic transition, in agreement 
with the numerical simulations of Selinger and Ratna Fig. H] shows, for zero applied 
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stress {a = 0), the curve in the 7-c plane above and to the right of which there is no 
first-order transition. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



c 

FIG. 4. The curve for a = above and to the right of which there is no first-order transition. The 
vertical dotted line corresponds to c = f /3. 

C. Results for an isotropic disorder (c = f/3, 7 > ) 

We recall that an isotropic distribution of random fields c = 1/3 represents samples 
formed in the isotropic state, with the random stresses coming from the cross-linked network 
without preferred direction. According to Fig. HI in this case a first-order transition occurs 
in zero applied stress (cr = 0) for any value of the random-field strength 7. To illustrate this 
fact, we plot in Fig. [5] the distortion factor A as a function of temperature in the absence of 
external stress {a = 0) for several values of the disorder parameter 7. The nematic phase 
decreases with increasing 7, but the first-order transition persists showing no evidence of 
the experimentally observed smoothed out nematic-isotropic transition in samples formed 
in the nematic state. 

In the presence of the applied stress (cr > 0), the first-order transition is smoothed out by 
a sufficiently large temperature or disorder strength. This can be seen in the stress-strain 
curves, where the strain e is related to the distortion by the equation A = 1 -|- e. In Fig. [6^ 
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FIG. 5. Distortion factor A as a function of temperature for a = 0, c = 1/3, and 7 between and 1. 

we plot isotherms for 7 = 0.6 and in Fig. |6]d we plot iso-7 curves for T = 0.8. It is clear 
from these figures that the temperature and the disorder strength have similar effect on the 
system. For sufficiently low temperature (T < Tc) or disorder strength (7 < 7c), the stress 
grows monotonically with strain in the isotropic phase up to the first-order transition to the 
nematic phase. At the transition the two phases coexist and the strain is independent of 
stress. The stress-strain "plateau" is then followed by a stress growth at larger strains. This 
behavior agrees with the experimental scenario of a typical transition between polydomains 
and a monodomain in NEs. As the temperature or the disorder strength increases, the 
coexistence curve shrinks until the critical point is reached at T = or 7 = 7^. The whole 
phase diagram resembles a typical P-V diagram of simple fluids. 
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V. DISCUSSIONS 



Recent experiments have shown that NEs crosshnked in the isotropic state display a 
well-defined plateau at the stress-strain curve, for substantially lower critical stresses, in 
comparison with NEs crosslinked in the nematic state [2^ . Optical microscopy observations 
suggest this behavior be attributed to larger memory effects for NEs formed in the nematic 
state. We may then expect that 7 is not independent, but should increase with c, according 
to the language of our model. Now it is not difficult to find appropriate values of 7 and 
c satisfying this restriction, and in general agreement with the experimental results. We 
show in Figure [7] two stress-strain curves for T = 0.95. For curve (a), we consider samples 
crosslinked in the isotropic state (c = 1/3), which implies no memory effect (7 = 0). In this 
case, the stress-strain curve presents a well-defined plateau for low critical stress. For curve 
(b), we consider samples crosslinked in the nematic state (c = 0.35), which should imply 
strong memory effects (7 = 0.45). Note that the stress-strain curve characterizes a broad 
polydomain-monodomain transition for hig her critical stress, in good agreement with the 
experimental results by Urayama et al. 



0.16 




e 



FIG. 7. Stress-strain curves for T = 0.95. (a) 7 = 0.45 and c = 0.35. (b) 7 = and c = 1/3. 

In conclusion, we have introduced a simple mean-field lattice model to describe the be- 
havior of nematic elastomers. This model combines the Maier-Saupe-Zwanzig theory of 



liquid crystals [15|, ll8|, ll9| and the lattice version, due to Selinger and Ratna [5|], of the 
Warner- Terentjev theory of elasticity [l|. We performed detailed calculations for a large 
range of parameters, with the inclusion of the effects of a quenched distribution of random 
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fields. A stress-strain coexistence curve may be obtained for systems cooled below a freezing 
temperature, which is analogous to the P-V diagram of a simple fluid, with the disorder 
strength playing the role of temperature. Below a critical stress, the characteristic tie lines 
resemble the experimental stress-strain plateau, and may be interpreted as signatures of a 
polydomain-monodomain transition. In the monodomain case, wc show that random-field 
disorder may soften the first-order transition between nematic and isotropic phases, provided 
the samples are formed in the nematic state. Beyond general agreement with some previous 
findings, we hope our results may motivate further experimental work on the stress-strain 
coexistence curve of nematic elastomers. 
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